Long-range interacting rotators: connection with the mean— field approximation 
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We analyze the equilibrium properties of a chain of ferromagnetically coupled rotators which 
interact through a force that decays as r~ a where r is the interparticle distance and a > 0. Our 
model contains as particular cases the mean field limit (a = 0) and the first-neighbor model (a — > 
oo). By integrating the equations of motion we obtain the microcanonical time averages of both the 
magnetization and the kinetic energy. Concerning the long-range order, we detect three different 
regimes at low energies, depending on whether a belongs to the intervals [0,1), (1,2) or (2, oo). 
Moreover, for < a < 1, the microcanonical averages agree, after a simple scaling, with those 
obtained in the canonical ensemble for the mean-field XY model. This correspondence offers a 
mathematically tractable and computationally economic way of dealing with systems governed by 
slowly decaying long-range interactions. 



One of most important questions in statistical mechan- 
ics refers to the connection between dynamics and ther- 
modynamics: To what extent a suitable ensemble average 
allows to predict the time average of a physical observable 
performed by our instruments in the laboratory? Or, in 
other words, which are the mechanical specifications of 
those systems to which the results of statistical physics 
can be applied? Within this context, while ergodicity and 
mixing have been analyzed intensively in the literature, 
there is another important point which has not deserved 
the same degree of attention, that is the possibility of 
defining a thermodynamically suitable energy function. 
In fact, for systems governed by sufficiently long-range 
interactions decaying as r~ a with the interparticle dis- 
tance r, there results a non-extensive Hamiltonian, i.e., 
the energy per particle diverges in the thermodynamics 
limit N — > oo Jj],||. Gravitational (a — 1, d = 3) and 
monopole-dipole (a = 2, d = 3) interactions are only two 
well known instances among many others. Furthermore, 
such forces are particularly interesting since they can lead 
to equilibrium behaviors different from those observed in 
short-range systems and even give place to phase transi- 
tions otherwise absent, even in the d = 1 case. 

Our aim here is to investigate how to deal with sys- 
tems governed by long-range interactions by analyzing a 
simple but rich prototype with adjustable a. The main 
goal of this letter is to show that the mean-field limit 
(a = 0) is able to describe the thermodynamics in the 
whole range < a < 1. The model consists in a one 



dimensional chain of N interacting rotators with peri- 
odic boundary conditions. Each rotator moves on the 
unit circle and therefore it is fully described by the an- 
gle — 7r < Oi < 7r and its conjugate momentum pi (with 
i = 1, . . . , N). The dynamics of the chain is governed by 
the following Hamiltonian 
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K + U, (1) 



where, without loss of generality, we have chosen unitary 
moments of inertia for all the particles. Here measures 
the minimal distance between rotators i and j along the 
chain. The Hamiltonian ([j]) describes a classical inertial 
XY ferromagnet. K and U denote the kinetic and poten- 
tial energies respectively. The equations of motion ruling 
this dynamical system are: 



t = Pi 

p i = -'£wx(0 i -0 j )/r% 



(2) 
(3) 



We associate to each particle a spin vector 

rrij = (cos 6»i, sin 6>i) (4) 
and define the total magnetization of the system as: 

M = ^E m - ( 5 ) 
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The long-time behavior of M determines whether the 
system orders (M / 0) or not (M = 0). 

We also introduce a time dependent temperature T(t) 
as T(t) = (2/JV) < K > (t), where < . . . > denotes a 
time average performed over the time interval (0,t). By 
calculating the long-time behavior of T as a function of 
the total energy of the system one gets the caloric curve 
T vs. E/N from which the specific heat function is ex- 
tracted. 

Note that the a — > oo limit yields the first -neighbor 
case while a = represents the mean-field version. The 
later case has a correspondence with the model known 
in the literature as hamiltonian mean-field XY (HMF), 
provided the potential energy (thus the strength of the 
interactions) is scaled by the number of particles N, 
i.e., in the HMF model U = ^ Ei/jl 1 ~ cos (°i ~ 
The HMF model has received special attention during 
the last years and very interesting results have been 
obtained both for its equilibrium and non-equilibrium 
properties | §. This model can be solved analytically 
within the canonical ensemble formalism ||. It is found 
that there is a critical specific energy at which the sys- 
tem suffers a phase transition separating a paramag- 
netic high energy phase (with M = 0) from a ferro- 
magnetic low energy one (with M/0). More precisely, 
E/N = T/2 + (I - M 2 )/2, with |M| = yT where y is the 
value maximizing — y 2 T/2 + In Io(y), being Iq the modi- 
fied Bessel function of order 0. For E/N > 3/4 (T > 1/2) 
the only solution is y = while below the critical energy 
the null solution becomes unstable and a new stable solu- 
tion appears giving place to a second order phase transi- 
tion. On the opposite limit of first-neighbor interactions 
(a — > oo), when U = ~ cos(6* i+ i — 0,)], one does 

not observe an order-disorder thermodynamical transi- 
tion @. In this case, E/N = T/2 + 1 - h(l/T)/I (l/T), 
where I n stands for the modified Bessel function of order 
n. The dynamics generated by the Hamiltonian (^) has 
been analyzed before [Q. In the limit N — » oo and for 
energies above a critical value, the maximal Lyapunov 
exponent was shown to vanish for < a < 1, while it 
tends to a finite value otherwise. 

For arbitrary values of a we have integrated numeri- 
cally the set of equations of motion (2) and (3) using a 
fourth order symplectic method [H with a fixed time step 
selected so as to keep the energy constant within an error 
AE/E of order 1CP 4 . Initial configurations (t = 0) were 
chosen as follows: all the angles were set to zero and the 
momenta were chosen at random from the uniform dis- 
tribution with zero mean value. Next, all the momenta 
were scaled in order to attain the desired total energy E. 
We varied the size of the systems from N — 100 to 1600 
in order to analyze finite size effects Q. Once elapsed a 
transient (that depends both on the size of the system 
and on the total energy E), we computed the time av- 
erages of both the magnetization and the kinetic energy 
looking for their asymptotic behavior. 



We start considering the region < a < 1. Systems 
characterized by values of a within that interval do not 
have a well defined Hamiltonian in the sense that, in 
the thermodynamics limit N — > oo, the potential energy 
diverges. One of the goals of this work is precisely to 
establish a connection between our microcanonical tem- 
poral averages obtained for a inside that range and those 
obtained theoretically by averaging in the canonical en- 
semble the HMF model ||. 

Fig. la displays the modulus of the magnetization |M| 
as a function of the energy per particle E/N for a = 0.5 
and different system sizes. We observe that the system 
undergoes a phase transition from a ferromagnetic low 
energy phase with M ^ to a paramagnetic high energy 
phase where M = 0. As expected, the critical specific 
energy diverges as N increases, due to the lack of exten- 
sivity of the system. In Fig. lb we present the same 
results but now with the specific energy scaled by a fac- 
tor N that depends both on the system size N and on a. 
To understand the origin of this scaling, let us stress that 
the lack of extensivity in the range < a < 1 emerges as 
a consequence of the divergence of the specific potential 
energy upper bound in (Jl|) . N is nothing else than twice 
the value of this upper bound: 



N/2 



N = 2V- 



fN/2 
2 drr~ 

Jl/2 



TV 1 "" - 1 

= 2" — , (6) 

1 — a 



For N — > oo, N behaves as 

( 2 a jz^N 1 -" for 

N(a) ~ \ 2 In TV for 

I 9(a) ^ for 



< a < 1 
a = 1 
a > 1 



(7) 



where O(a) is a function of a which, for a — > oo goes to 
2 (a — 1). Note that all the curves for different iV collapse 
into a unique one, despite small discrepancies around the 
critical value. Similar TV-scaling collapse was already ob- 
served for magnetic systems 0,|ll[ and also for systems 
governed by interactions of the Lennard- Jones type p2[ . 
In all these cases, although the a-dependent prefactor 
in (^) is model dependent, the behavior of N with N is 
invariant. 

Fig. lc plots T/N vs. the scaled specific energy 
E/(NN) (caloric curve) also for a = 0.5 and the same 
system sizes considered in Figs, la and lb. Here again 
the iV-scaling leads to data collapse. It is worth here to 
stress that around the critical energy this plot depends 
strongly on the equilibration transient, the size of the 
system and the initial configuration adopted. Our re- 
sults within this parameter range seem to indicate the 
existence of a first order phase transition, with the high 
energy phase coexisting with the ordered one. An analo- 
gous behavior, already reported for the HMF model [0, 
where the transition is second order, was believed to be 
a purely microcanonical result reflecting the existence of 



2 



long living quasistationary nonequilibrium states whose 
lifetimes increase with N. This seems to be also the case 
for any < a < 1 since the discrepancy around the 
transition for fixed energy and fixed size is attenuated by 
averaging over larger time intervals. 

Fig. 2 exhibits |M| vs. Ej (NN) for a = 0.25, 0.5 and 
0.75, all for N = 400. Observe that the curves agree with 
a unique one! The same collapse is also detected for the 
caloric curve (not shown). Let us recall the main motiva- 
tion of this letter, namely the possible relation between 
statistical and temporal averages for systems for which 
one cannot a priori define an extensive energy, such the 
case we are considering now. In Fig. 2 we have also 
included the plot (solid line) of the theoretical predic- 
tions of the equilibrium values obtained, by means of 
the canonical ensemble formalism, in the HMF version 
(a = with TV-scaled potential energy) ||. What we 
now observe is that not only different size and different 
a curves, with < a < 1, collapse into a unique one, but 
they collapse precisely to the mathematically tractable 
extensive mean-field limit! A similar collapse had al- 
ready been found for the Ising ferromagnet submitted to 
a Monte Carlo process @ but this is, as much as we 
know, the first time this effect can be confirmed in a 
conservative model with deterministic dynamics. Since 
slowly decaying long-range interacting systems are ubiq- 
uitous in nature, our results, if valid for all those systems 
such that a < d, reveal a simple way of calculating quan- 
tities at the equilibrium. 

Next we briefly describe what happens for 1 < a < 2. 
There exists still an order-disorder transition but with 
a sensitive dependence on the value of a, as can be ob- 
served in Figs. 3a and 3b for the particular case a = 1.5. 
Contrary to the previous case where the critical energy 
(above which the magnetization falls down to zero as N 
increases) remains independent on a although slightly 
smaller than the theoretical prediction, here as the value 
of a approaches 2, the critical energy decreases until a 
finite value which can be non-zero |13|| . 

Finally, we analyze the region a > 2. Here, as ex- 
pected, the system behaves similarly as in the first- 
neighbor limit, i.e., it does never order since for any finite 
energy the magnetization goes down to zero as N — > oo. 
In Fig. 4 we plot T/N vs. E/(NN) obtained numeri- 
cally (symbols) for a = 2.5 and a = 5.0, together with 
the theoretical values for both the limit a — > oo [§ (solid 
line) and the HMF model || (dashed line). We see how 
the numerical results approach those of the a — > oo limit. 

Summarizing, we have studied the equilibrium behav- 
ior of a one-dimensional conservative system of interact- 
ing particles as a function of the range of the interactions 
a. By integrating numerically the equations of motion 
we have found three different classes of systems. For 
< a < 1 the systems undergo a second order phase 
transition and the measured quantities (e.g., the critical 
energy) when suitably scaled, do not depend on the value 



of a. For 1 < a < 2 the systems undergo also a sec- 
ond order phase transition but, unlike the previous case, 
the critical energy and the magnetization curve depend 
sensitively on the range of the interactions a. Finally, 
for a > 2 the systems adopt the first-neighbor behavior 
where there is no order at finite temperature for N — > oo. 

It is also worth noting that an the iV-scaling performed 
over the results produced by Hamiltonian leads to 
the same results as those produced by an artificial "ex- 
tensive" Hamiltonian constructed with iV-scaled poten- 
tial energies generalizing the HMF approximation. For 
a correspondence between both treatments when looking 
at dynamical aspects, time should be A^ 1 / 2 -scaled. 

Probably the non-trivial data collapse here reported 
for < a < 1, can be extended to < a < d for an 
arbitrary dimension d, being (from the generalization of 
(|) to d dimesions) N ~ T^/d{N 1 ~ a/d - 1). If true, our 
findings reveal a simple way of calculating macroscopic 
quantities at the equilibrium when long-range interac- 
tions are involved. 
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CAPTIONS FOR FIGURES 

Figure l:(a) Modulus of the magnetization, |M|, as a 
function of the specific energy E/N; (b) |M| as a function 
of the scaled specific energy E/(NN); (c) scaled temper- 
ature T/N as a function of E/(NN). The symbols cor- 
respond to numerical simulations in the microcanonical 
ensemble for a — 0.5 and different system sizes indicated 
on the figure. Each symbol corresponds to an average of 
different initial conditions (typically 10). The solid lines 
correspond to the theoretical equilibrium predictions for 
the HMF model (for which N = 1, since the HMF Hamil- 
tonian is already TV-scaled). 

Figure 2: |M| as a function of E/{NN). The symbols 
correspond to numerical simulations in the microcanoni- 
cal ensemble for a = 0.25, 0.5, and 0.75, all for N = 400. 



The solid line corresponds to the theoretical results for 
the HMF model. 

Figure 3: (a) |M| as a function of E/(NN); (b) scaled 
temperature T/N as a function of E/(NN). The sym- 
bols correspond to numerical simulations in the micro- 
canonical ensemble for a — 1.5 and different system sizes. 
For comparison, the dotted line corresponds to the theo- 
retical equilibrium predictions of the HMF model. 

Figure 4: Scaled temperature T/N as a function of 
E/(NN). The symbols correspond to numerical simula- 
tions in the microcanonical ensemble for a = 2.5 (gray), 
5.0 (white) and different system sizes indicated on the 
figure. The dashed and solid lines correspond to the the- 
oretical results for the HMF (N = 1) and the a — > oc 
(N = 2) limits, respectively. 
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